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Abstract. The formation of a double white dwarf binary likely involves a common 
envelope (CE) event between a red giant and a white dwarf (WD) during the most 
recent episode of Roche lobe overflow mass transfer. We study the role of recom¬ 
bination energy with hydrodynamic simulations of such stellar interactions. We find 
that the recombination energy helps to expel the common envelope entirely, while if 
recombination energy is not taken into account, a significant fraction of the common 
envelope remains bound. We apply our numerical methods to constrain the progenitor 
system for WD 1 IOIh- 364 - a double WD binary that has well-measured mass ratio of 
q - 0.87 + 0.03 and an orbital period of 0.145 days. Our best-fit progenitor for the 
pre-common envelope donor is a 1.5 Mq red giant. 


1. Introduction 

The formation channel of a double white dwarf (DWD) binary is still not fully certain, 
although it is believed that the formation of a compact binary system composed of two 
white dwarfs (WDs) includes a common envelope (CE) event, at least during the last 
episode of mass exchange between the first-formed WD and a low-mass giant. Low- 
mass giants have a well-defined relation between their core masses and radii. This 
makes the state of the progenitor binary system at the onset of the CE event leading to a 
DWD the best theoretically understood from among the progenitors of all known types 
of post-CE systems. However, previous attempts to model a CE event between a low- 
mass giant and a WD did not succeed to eject the entire CE during three-dimensional 
hydrodynamical simulations - a significant fraction of the expanded envelope remained 
bound to the formed binary, with almost no energy transfer taking place from the binary 
orbit to that circumbinary envelope (Passy et al. 2012; Ricker & Taam 2012). It has been 
proposed a while ago that the recombination energy of H and He should play a role in 
unbinding the donor’s envelope (Lucy 1967; Roxburgh 1967; Paczyhski & Ziolkowski 
1968; Han et al. 1994, 2002); however, to date, this energy has not yet been taken into 
account in the modeling of a CE event. 
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2 . Numerical method 

We use a 3D Smoothed Particle Hydrodynamics (SPH) code, StarSmasher (Gaburov 
et al. 2010; Lombardi et al. 2011); for details on modeling of CE events see Nandez 
et al. (2014). This code by default uses a standard analytical SPH equations of state 
(EOS) that takes into account radiation pressure and ideal gas pressure, but not the 
changes in gas’ ionization states. We added a possibility to run the code using a tab¬ 
ulated EOS adapted from MESA (Paxton et al. 2011), where recombination energy is 
taken into account. To model donor stars, we first evolve one-dimensional low-mass 
stars using EV/stars (recent updates described in Glebbeek et al. 2008) and then relax 
them in StarSmasher. The comparison of the initial relaxed stars with two different 
EOSs have shown that SPH profiles in bofh cases mafch very well fhe pressure, densify, 
and gravitational potential with the stellar profiles. However, fhe SPH specific internal 
energy profile is matched only in case when the tabulated EOS is used (see Eigure 1). 



Figure 1. The specific internal energy m as a function of the radius coordinate, for 
the I.SMq giant with 0.31 Mq core. The solid line corresponds to one-dimensional 
stellar profile, the gray triangles correspond to the relaxed star with the tabulated 
EOS, and the black boxes are for the standard EOS. 


3. Initial conditions 

One of the best measured DWD systems is WD 1101-1-364. This system has an orbital 
period of 0.145 d and the mass ratio between the younger and older WDs q = mi/ni 2 = 
0.87 ± 0.03, where nii = 0.3 IM© and m 2 - 0.36Mq (Marsh 1995). Accordingly, we 
adopted for our CE simulations that the WD companion has a mass of O.36M0. Eor a 
red giant companion, we considered low-mass stars with masses 1, 1.2, 1.4, 1.5, 1.7 and 
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1.8 Mq. We have evolved those stars until the moment when the core has a He mass of 
0.31 - O.32M0. Figure 2 shows how the stellar radius changes when the He core mass 
grows. We adopted that the radius that a star has when its core mass is 0.31 - 032Mq 
is equal to the Roche lobe radius defined as in Eggleton (1983). This provided us with 
the initial orbital separation of the progenitor binary at the start of our CE simulations. 
Because the radii of the giants are a strong function of the core mass and not the total 
mass, the initial orbital separations do not vary strongly among the cases with giants of 
different masses. 



Figure 2. Stellar radii as a function of the red giant’s He core mass. The stellar 
tracks correspond to 1.0, 1.2, 1.4, 1.5, 1.7 and 1.8 Mq. The gray vertical lines in¬ 
dicate the region of possible initial radii of the red giants, to be consistent with the 
observation of WD 1101-r364. 


4. Results 

We compare the two EOSs for the example of an initial binary consisting of ESM© giant 
star (modeled with 99955 particles) and a O.36M0 WD (modeled as a point-mass). Eor 
each particle, we find its total energy as the sum of its kinetic, potential and internal 
energies (internal energy here includes recombination energy, if still stored). We then 
define: 

• Ejecta - fhe unbound material where each particle has its total energy positive; 

• Binary - all the particles that have their total energy negative, and are located 
inside of one of the two Roche lobes; 
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• Circumbinary - all particles which are bound to a formed binary (the total energy 
of each particles is negative), and are located outside of the either of the Roche 
lobes. 


The simulations were stopped after 6000 orbital periods counted from the compact 
binary formation. See the Table 1 for the masses and total energies at the end of the 
simulations. In our simulations, just as in the previous studies (Passy et al. 2012; Ricker 
& Taam 2012), the run with the standard EOS resulted in more than a half of red giant 
envelope to remain bound to the system, even though the value of its total energy is 
very small compared to the orbital energy of the binary. At the same time, the run that 
used the tabulated EOS resulted in the whole envelope being expelled. 

The energy formalism that described a CE event simplistically (for an overview, 
see Ivanova et al. 2013) predicts that the more massive the donor is, the tighter the orbit 
gets. Indeed, we find the same trend with our simulations - see Eigure 3, where we 
show how the final orbital period changes with the initial giant’s mass, only for the runs 
that were performed using the tabulated EOS. 

Our best match with the observed orbital period P - 0.145d (Marsh 1995) is 
obtained for a donor of 1.5 M©. 



Figure 3. Final orbital periods, as a function of the pre-CE giant mass. The dotted 
and dashed curves represent, respectively, simple linear and quadratic fits to our 
simulation result. The solid line represents the observed period of WD 1101+364. 
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Table 1. CE event involving 1.5Mq giant and 0.36Mq WD. 


EOS 

^^unb 

«Uir 

^^bin 

poo 

tot.unb 

tot,cir 

poo 

tot,bin 

standard EOS 

0.521 

0.658 

0.679 

2.150 

-0.402 

-24.011 

tabulated EOS 

1.178 

0.000 

0.680 

6.910 

0.000 

-25.508 


niunb is the mass of the ejecta, nicir is the mass of circumbinary material, and mbin is the 
mass of the formed binary, all in Mq. ^Mcir’ ^totbin energy of the 

ejecta, circumbinary material, and the binary, respectively, all in units of 10"^^ erg. 

5. Conclusion 

With our simulations, we have confirmed that the recombination energy plays an im¬ 
portant role in the ejection of a CE, where it helps to unbind the entire envelope. 

The observations for WD 1101-1-364 are well reproduced with a progenitor binary 
composed of a l.SM© red giant and a O.ShM© WD companion. In this case, the modeled 
post-CE binary has a final orbital period of Porb - 0.135 d, the newly formed WD 
has a mass of O.319M0 and an older WD has a mass 0.361Mo; binary’s mass ratio is 
q = 0.884, in agreement with the observational error. 
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